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INDUCED VELOCITIES NEAR A LIFTING ROTOR 
WITH NONUNIFORM DISK LOADING 1 


By Harry H. Heyson and S. Katzoff 


SUMMARY 

A theoretical and experimental investigation of the induced 
velocities near a lifting rotor has been conducted. The results 
are presented in three parts. 

Part I presents a method for computing , from available calcu- 
lations for uniform disk loading , the effect of nonuniform 
circularly symmetrical disk loading on the normal component 
of induced velocity in the vicinity of a lifting rotor. Charts of 
the normal component of induced velocity are given for the longi- 
tudinal plane of symmetry and for the major axes of rotors with 
two different , nonuniform circularly symmetrical disk-load 
distributions. It is shown, that the normal component of in- 
duced velocity must be zero at the center of any practical rotor . 
A comparison of the results of this report with those for a uni- 
form disk loading shows that nonuniform disk loading has a 
powerful effect on the induced velocity distribution and that it 
must be taken into account in estimating the effect of the rotor 
on most components of an aircraft. 

Part II develops certain symmetry relations for the induced 
velocities in the plane of a uniformly loaded rotor and also de- 
velops relations between the radial load distribution of the rotor 
and the radial variation of induced velocities in the wake. 

Part III presents the results of an investigation in the Langley 
full-scale tunnel of the induced flow near a lifting rotor. Meas- 
urements of stream angles and velocities were made in several 
transverse planes along and behind the rotor in four different 
conditions representative of the cruising and high-speed ranges 
of flight. These measurements indicate that available theory 
may be used to calculate with reasonable accuracy the induced 
flow over the forward three-quarters of the disk for these flight 
conditions provided that a realistic nonuniform rotor disk-load 
distribution is assumed. Rearward of the three-quarter- 
diameter point, calculations of the induced velocity are increas- 
ingly inaccurate due to the rolling up of the trailing -vortex 
system. Farther rearward , well behind the rotor , the flow may 
be represented more accurately by the flow behind a uniformly 
loaded wing. 

INTRODUCTION 

A knowledge of the flow induced by a lifting rotor is neces- 
sary in any detailed study of the characteristics of a heli- 
copter or convertiplane. The induced flow is one of the most 


important factors involved in fuselage downloads or in the 
performance and loading of the wings of convertiplanes. 
Large effects on stability may also be experienced as a result 
of the induced flow acting on the tail surfaces. The rotor- 
induced flow actually involves a problem of unsteady aero- 
dynamics, however, a knowledge of the time-averaged ve- 
locities should be sufficient for study of the problems men- 
tioned. Experimental evidence (ref. 1) indicates that the 
time-averaged induced velocities may also be sufficient to 
estimate the interference between the rotors of multiple- 
rotor configurations. However, the differences between the 
time-averaged and the instantaneous velocities (ref. 2) is 
such that a more rigorous approach is probably required for 
accurate analyses (for example, blade vibrations) of the 
rotor system. 

Even the time-averaged velocities have proved difficult to 
compute, but several simplified analyses are available (for 
example, refs. 3 to 6). These analyses are limited by the 
restricted number of locations for which the calculations 
have been made and also by the fact that the loading is 
usually assumed to be uniform. The most complete treat- 
ment is that of reference 6, which computes the normal com- 
ponent of induced velocity along the rotor axes and through- 
out the longitudinal plane of symmetry. Despite the as- 
sumption of uniform loading and the restriction to a single 
plane, it was hoped that the calculations would prove ade- 
quate to explain the effect of the rotor on short-span tail 
surfaces. Such, however, is not the case, for the results of 
wind-tunnel tests of the present investigation indicate that 
the induced velocity at typical tail locations may differ 
from that of reference 6 by as much as 1 y 2 times the mean 
induced velocity at the rotor. One principal cause of this 
difference is the fact that the disk loading on an actual rotor 
is far from uniform as is assumed in reference 6. 

In view of the difficulty of calculating the induced veloci- 
ties, it would be expected that comprehensive measurements 
of these velocities would be available. However, the actual 
direct measurement of these velocities in flight or in small 
wind tunnels has proved to be exceedingly difficult. Several 
previous investigations have provided some qualitative in- 
formation by the use of smoke, but quantitative measure- 
ments are relatively scarce. (A comprehensive review and 


1 Supersedes NACA Technical Note 3690 by Harry H. Heyson and S. Katzoff, 1956, and 3691 by Harry H. Heyson, 1956. 
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bibliography on induced-flow measurements is given in ref. 
7.) Such measurements as are available are generally inade- 
quate in scope to provide designers with sufficient experi- 
mental information or to define the general nature of the flow. 

The present investigation attacks the problem both the- 
oretically and experimentally. The report is divided into 
three parts as follows: Part I uses the results of reference 6 
to build up by superposition the induced flow near a rotor 
having any arbitrary axisymmetric nonuniform distribution 
of disk load. For two rotors having different realistic disk- 
load distributions, charts of the normal component of induced 
velocity are given. Part II contains analytical treatments of 
some interesting characteristics of the downwash field that 
were observed during the course of the calculations for part 
I. More specifically, this part develops certain symmetry 
relations for the induced velocities in the plane of the uni- 
formly loaded rotor as well as relations between the radial 
load distribution of the rotor and the radial variation of 
induced velocities in the wake. Part 111 presents the results 
of an experimental investigation conducted in the Langley 
full-scale tunnel. Stream-angle and dynamic-pressure meas- 
urements were made in several transverse planes near a 
rotor. The simulated flight conditions were chosen so as to 
represent a helicopter in level flight at normal cruising 
speeds. The tip-speed ratios range between 0.095 and 0.232; 
the calculated skew angles, between 75.0° and 85.8°. Essen- 
tially all of the basic data are presented in order to provide 
the maximum possible amount of experimental information 
for design purposes. At those locations for which theoretical 
calculations are available, the data are compared with theory 
in order to indicate the extent to which the theory is ade- 
quate and also to determine the general nature of the flow. 

SYMBOLS 


The following symbols are typical of usual rotary -wing 
terminology and are used throughout parts I and 111 of 
this report. The usual complex-variable symbols are used 
in part 11 since they are more convenient for that analysis. 
The symbols for part II are defined separately in that portion 
of the report. 

A area 


b 

r L 

Cr 

I 


i 


L 
< 2 
<l.o 


number of blades 
rotor lift coefficient, 


L 



T 


rotor thrust coefficient, 

p (ft/r) 7r/c z 

equivalent flat-plate area representing parasite drag, 
based on unit drag coefficient, 

Helicopter parasite drag 

pK! 

2 


sq ft 


complement of skew angle (used in ref. 5 as angle 
of attack of rotor tip-path plane), deg 
lift, lb 

dynamic pressure, lb/sq ft 
free-stream dynamic pressure, lb/sq ft 


r blade-element radius, ft 

R rotor radius, ft 

R v radius of vortex cylinder, ft 

s semispan of wing, ft 

T rotor thrust, lb 

v normal component of induced velocity, positive 

downward, ft/sec 

Vq average or momentum theory value of the normal 
component of induced velocity, positive down- 
ward, ft/sec 

V forward speed of rotor, ft/sec 

X distance parallel to longitudinal rotor tip-path-plane 
axis, measured positive rearward from center of 
rotor, ft (fig. 1) 

Y distance parallel to lateral rotor tip-path-plane axis, 

measured positive on advancing side of disk from 
center of rotor, ft (fig. 1) 

Z distance parallel to axis perpendicular to rotor tip- 

path plane, measured positive above rotor from 

center of rotor, ft (fig. 1) 

a rotor control-axis angle of attack, radians 

7 vorticity, radians/sec 

T circulation, ft 2 /sec 

r' circulation at tip of blade, ft 2 /sec 

x - r , 4 . V sin ol—v o 

X inflow ratio, 

. . V cos ol 

i u tip-speed ratio, — ^ — 

p mass density of air, slugs/cu ft 

X rotor wake skew angle, angle between Z-axis of 

tip-path plane and axis of skewed wake, positive 
rearward from Z-axis, deg (fig. 1) 

£2 rotor angular velocity, radians/sec 

Note that the definitions of X, Y , and Z, together with 
the direction of rotation of the rotor in part III of this in- 
vestigation, result in a left-handed axis system. 



Figure 1. — Coordinate system of rotor and wake. Arrows denote 
positive direction. 
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I. DEVELOPMENT OF INDUCED VELOCITY FIELDS BY SUPERPOSITION 

By Harry H. Heyson 


THEORY 

THE WAKE CONSIDERED AS A SYSTEM OF COAXIAL VORTEX CYLINDERS 

For the present simplified theory, the forward speed of 
the rotor is assumed negligible relative to the rotational 
speed of the blades. Then, at any radial location, the 
local velocity of a blade element is simply fir, so that the 
lift per unit radial distance can be written 

~=bp9.rT ( 1 ) 


where T, in general, is a function of r. The local disk 
loading (or annulus loading) is 


Local disk loading 


bpttrT dr 6pfiT 

2i rr dr 2 t 


(2) 


In general, as illustrated by equations (1) and (2) } the 


local disk loading varies as - times the blade loading. Thus, 

a uniform disk loading implies a triangular blade loading, 
and a triangular disk loading implies a parabolic blade 
loading. 

Equation (2) also shows that for a uniform disk loading 
the circulation is constant along the blade. For this case, 
vorticity is accordingly shed only at the blade tips, and the 
wake consists of a single vortex cylinder (fig. 2 (a)). In 
references 3 and 6, where the lifting rotor is represented by 
a uniformly loaded disk, the downwash field is correspond- 
ingly calculated as the downwash field of this single vortex 
cylinder. 

If, however, the disk loading is not uniform, the circulation 
varies along the length of the blade and vorticity must be 
shed all along the blade instead of only at the tip. Con- 
sequently, a rotor having nonuniform disk loading must be 
represented by a distribution of shed vortex cylinders. 

In the present calculations, these cylinders are all assumed 
parallel and concentric (see fig. 2 (b)) and are inclined at a 
skew angle defined by the forward velocity of the rotor 
and the momentum-theory value of the rotor induced 
velocity. Such an assumption presents obvious incon- 
sistencies, because such an idealized geometry could not 
be produced or maintained in the nonuniform induced-flow 
field that it itself creates. In spite of such inconsistencies, 
however, the arrangement of concentric cylinders is probably 
a better basis for computations than a single vortex cylinder 
trailing from the blade tips. 

It is possible to consider that the solution of reference 6 
represents the flow generated by any one of these vortex 
cylinders. The radius used in that report should now be 
considered to be the radius R v of the individual vortex 
cylinder. The contribution of any single vortex cylinder 


to the total induced flow at a point 



may be read 


directly from the charts of reference 6 at the point 

[{XIR)(R/R,), (Y/R)(R/R V ), (Z/R)(R/R v )l 




(a) Uniform load. 

(b) Nonuniform axially symmetric load. 


Figure 2. — Assumed vortex pattern of rotor wake. 

In practice, any arbitrary circularly symmetrical rotor 
disk-load distribution may be approximated by the sum of 
the loads carried by a finite number of superposed con- 
centric vortex cylinders of suitable strength and dimensions. 
Then the induced velocity field of the nonuniformlv loaded 
rotor will be the sum of the superposed induced velocity 
fields of the vortex cylinders comprising its wake. 

Some difficulty may be anticipated because of the dis- 
continuities in the flow at the edges of the vortex cylinders. 
In practice, however, only a small error occurs if a sufficient 
number of vortex cylinders are used, and the points of dis- 
continuity themselves are avoided. 

SAMPLE CALCULATION 

In order to demonstrate the method of superposition, the 

induced velocity at the point 0^=0 .2,^=0, ^=0.0 will be 

calculated for a rotor with a triangular disk loading and 
operating at a skew angle of 90°. 

As shown previously, the local disk loading is proportional 
to T; therefore, for uniform disk loading the circulation is 
constant along the blade. If the disk loading is triangular, 
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so that the local disk loading is 

Local disk loading= r' - 

— 7T 21 

and the average disk loading for the entire rotor is 

Average disk loading=— f r' ~ 2 wrdr 

7r/rJo 2 tt h 

b P nr' 

Sir 

Thus, for this case of the triangular disk loading, 

Local disk loading=| pX (Average disk loading) 


i 


i 


i 


Thus, the vortex field of a rotor with a triangular loading 
may be represented by one positive vortex cylinder with the 
same radius as the blade tips and a strength 1.5 times as 
great as a uniformly loaded rotor, plus ten equal negative 
vortex cylinders of smaller radii and a strength 0.15 times as 
great as a uniformly loaded rotor. The sum of the resulting 
disk loads is shown in sketch 1 and adequately represents a 
triangularly loaded rotor: 



Sketch 1. 

The contribution of each vortex cylinder may be found 
from figure 4 (g) of reference 6. 

The computation may be carried out in tabular form as 


follows: 

(§■) 

©(#.) 

©(£> 

(IX#.) 

( £ ) 

\Vo/ v 

1. 00 

0. 200 

0 

0. 100 

1. 15 ( 

1. 50) = 

1. 725 

. 95 

. 211 

0 

. 106 

1. 15 (- 

-. 15) = 

173 

. 85 

. 236 

0 

. 118 

1. 17 ( 

-. 15) = 

. 1 76 

. 75 

. 267 

0 

. 134 

1. 19 (- 

-. 15) = 

-. 179 

. 05 

. 308 

0 

. 154 

1. 20 (- 

15) = 

-. 180 

. 55 

. 364 

0 

. 182 

1. 21 (- 

-. 15) = 

. 182 

. 45 

. 445 

0 

. 223 

1.23 (- 

-• 15) = 

-. 185 

. 35 

. 571 

0 

. 286 

1. 24 (- 

-. 15) = 

186 

. 25 

. 800 

0 

. 400 

1. 25 ( 

15) = 

-. 188 

. 15 

1. 330 

0 

. 667 

. 95 ( 

-. 15) = 

. 113 

. 05 

4. 000 

0 

2. 000 

. 20 (• 

-. 15) = 

039 





v/v„= 

z(v/v o) v - 

= 0. 094 


It will be noted, in general, that the charts of reference 6 
must be extrapolated to obtain the contribution of some of 
the smaller cylinders. Also, if the maximum loading does 
not occur at the tips, more than one positive vortex cylinder 
will be required to adequately represent the flow. 


Sketch 2. 

INDUCED VELOCITY AT CENTER OF ROTOR 

Consider a rotor with a disk loading that is uniform 
everywhere except in a central “cut-out” portion where it 
carries no load. (See sketch 2.) The vortex field of this rotor 
may be represented by just two vortex cylinders of equal, 
but opposite, vorticity. Since the smaller cylinder has the 
same downwash field as the larger cylinder (except for 
dimensions), it follows that the induced velocity is zero at 
the center of the rotor disk. The same result, of course, 
will follow whenever any rotor has zero disk load at its 
center. Since all practical rotors have some cut-out near 
their centers (due to hub, etc.), they all must have zero 
induced velocity at their centers. 

This result is a special case of a more general theorem 
which is occasionally useful in comparing the results of 
theoretical calculations made for several different load 
distributions. Consider a rotor with a step loading com- 
posed of a number of uniform loads of different radii such 
as in sketch 3. The vortex field of this rotor may be repre- 
sented by several vortex cylinders of different strength. 
A uniformly loaded rotor produces an induced velocity 
of r 0 at its center. Therefore, the induced velocity at the 
center of the step-loaded rotor will be r 0l (due to uniform 
load 1) plus Vo 2 (due to uniform load 2) and so on. The 
final result is that the induced velocity at the center of the 
step-loaded rotor will then be the same as if the local load 
at the center existed over the entire disk rather than just 
at the center. 


- 


t 


Sketch 3. 
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A LIMITATION ON THE ANALYSIS 

It is important to note that the foregoing analysis is 
predicated upon a wake consisting of straight elliptic cylin- 
ders in which all of the vorticity is carried downward at some 
uniform velocity rather than the local velocities. The 
effect of this assumption is greatest in hovering for here 
the induced velocities are the solo velocities present. Since 
in this case the flow near the rotor center is carried away 
from the rotor at a much slower rate than the flow near the 
edges of the disk, the vorticity along the wake cylinders 
near the center must be correspondingly greater. In this 
extreme case, an appreciable error will be incurred in the 
computed induced velocities. However, when an apprecia- 
ble flight velocity is imposed on the rotor, as in rapid vertical 
climb or in forward flight, the vortices will be carried down- 
stream at a more uniform rate and the error will be small. 

Since the assumptions of the present analysis are grossly 
violated in static thrust, this case may be used to give an 
indication of the maximum possible error. From momen- 
tum theory, for any annulus 


dT=^'KVf>v 2 dr 


dT_ 

dA 


J_ dT 
2wr dr 


=2 pv 2 


which shows that in static thrust, the local induced velocity 
should be proportional to the square root of the local disk 
loading. In contrast, under the assumptions of the present 
report, the induced velocity has been shown to be directly 
proportional to the disk loading. The effect of this differ- 
ence on the induced velocities in the plane of a hovering 
triangularly loaded rotor is shown in figure 3. The momen- 
v / 3 /• 

turn theory gives — whereas the present vortex 

\ - it 

v 3 v 

theory gives — The maximum difference in induced 
Vo - it 

velocity occurs at the tip of the rotor and is 0.27bv 0 . Thus 
the results of the present vortex theory may be adequate 
to explain qualitatively certain trends even in hovering. 
For example, recent unpublished results from flow surveys 
below a hovering rotor indicate a region of up wash below 
the rotor center. Calculated results from the present 
theory for a triangularly loaded rotor also indicate the 
presence of an upwash in this region. 

In forward flight the assumptions of this report are much 
less severe. For example, reference 6 gives the mean in- 
duced velocity at the rotor as 

l QRC r 

v °~7 3 \ , 

(l-| m 2 ) a X 2 +m 2 

For tip-speed ratios greater than 0.10, X is small compared 
with n so that the induced velocity will be proportional to 
the loading as herein assumed. Thus the charts which will 
be presented in the following section should be sufficiently 
accurate for the rotor in forward flight. 



Figure 3. — Comparison of the results of momentum theory and present 
analysis for the induced velocities in the plane of a hovering, tri- 
angularly loaded rotor. 


CALCULATED RESULTS 

LONGITUDINAL PLANE OF SYMMETRY 

A number of nondimensional charts of the induced velocity 
have been computed. These charts not only illustrate the 
effect of a nonuniform disk loading, but they may also be 
used as working charts for the determination of the induced 
velocity at points in the field of practical rotors. The cal- 
culations were made for rotors having two different, circu- 
larly symmetrical, nonuniform disk loadings. These as- 
sumed disk loadings are shown in figure 4. The triangular 
loading was chosen since it is the simplest approximation to 
the actual rotor disk loading in all flight conditions. The 
variation of disk load labeled “typical load” was measured 
at a typical cruising condition during tests in the Langley 
full-scale tunnel of a rotor equipped to measure the rotor- 
blade pressure distribution. 

Figure 5 presents contour charts of the normal component 
of induced velocity in the longitudinal plane of symmetry 
for the rotor with a triangular disk loading for six different 
skew angles in the range from 0° to 90°. Figure 6 presents 
similar charts for a rotor having the typical measured disk- 
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Figure 4. — Disk loadings considered in calculations. 


load distribution for three different skew angles between 
63.43° and 90°, which bracket the skew angle at which this 
load distribution was measured. Figures 5 and 6 may also 
be used to find the induced velocity distributions for skew 
angles between 90° and 180° (the autorotation range) if 
the Z/R scale is multiplied by —1 and the chart for the sup- 
plement of the desired skew angle is read. 

LONGITUDINAL AXIS 

Figures 7 and 8 present the induced velocity distribution 
along the X-axis (or longitudinal axis) of the rotors with the 
triangular and typical measured load distributions, respec- 
tively. 

LATERAL AXIS 

Figures 9 and 10 present the induced velocity distribution 
along the 1 -axis (or lateral axis) of the same two rotors. 

Figures 7 to 10 may be used for skew angles from 90° 
to 180° merely by reading the curves for the supplement of 
the desired angle. 

Figure 9 indicates the order of accuracy of the calculations. 
Part II of this report will show that the variation of induced 
velocity on the lateral axis must be triangular if the disk 
loading is triangular. Actually, figure 9 indicates a slight 
upward curvature. The difference between these curves and 
straight lines is, however, small and could probably have 
been avoided by using a somewhat larger number of cylinders 
to represent the wake. 

OTHER REGIONS 

It would be interesting to carry out these calculations 
over the lateral planes where auxiliary devices such as wings, 
tail surfaces, or propellers would be located on flight machines. 
However, such calculations would be dependent upon a 
knowledge of the flow field of the uniformly loaded rotor and 
this flow field has not yet been computed. The flow in 
these regions can be calculated, if required, from considera- 
tions in reference 6 and the present report. 

DISCUSSION OF CALCULATED RESULTS 

LONGITUDINAL PLANE OF SYMMETRY 

Figure 1 1 presents a chart of the induced velocity in the 


longitudinal plane of symmetry of a uniformly loaded rotor 
at a skew angle of 63.43° (reproduced from ref. 6). 

The effect of nonuniform axisymmetric loading may be 
seen by comparing figure 11 with figures 5 (d) and 6 (a). 
Since both of the assumed nonuniform loads are zero at the 
center of the rotor, the induced velocity at the center is zero 
in both cases, as opposed to a value of v/v 0 =1.0 at the center 
of the uniformly loaded rotor. Both of the nonuniformly 
loaded rotors show an appreciable area of up wash just below 
and behind the center of rotation. They also show zero 
induced velocity further rearward on the center line of the 
wake. In this plane, the chart for the triangular load distri- 
bution shows a roughly triangular induced-velocity profile 
across the wake, and the chart for the typical measured load 
distribution shows a more rounded profile across the wake 
with a larger area of zero induced velocity reflecting the 
effect of the cut-out region in the vicinity of the hub. In 
contrast, the results for the uniformly loaded rotor case 
show the induced velocity building up continuously as the 
flow passes through the rotor and rearward in the wake until 
it reaches a uniform value of v/v 0 =2.0 near the trailing edge 
of the rotor. 

LONGITUDINAL AXIS 

Figure 12 (obtained from figs. 7 and 8 and ref. 6) shows the 
calculated effect of rotor disk-load distribution on the 
induced-velocity distribution along the X-axis (longitudinal 
axis) of the rotor. This figure clearly shows the large effect 
caused by the requirement that the induced velocity be zero 
at the center of the rotor for each of the nonuniformly 
loaded cases. Each of the nonuniformly loaded rotors has 
an upwash just behind its center, whereas the uniformly 
loaded rotor has downwash values of v/v 0 >1 in this region 
The differences in induced velocity between the nonuniformly 
loaded rotors are small along this axis, however, the rotor 
with the typical measured load distribution does have a 
somewhat greater upwash behind the center as a result of 
the larger cut-out region. 

LATERAL AXIS 

Figure 13 presents a similar comparison along the F-axis 
(lateral axis) of the rotors at a skew angle of 90°. This 
figure shows that the lateral center-line distribution is 
closely similar in shape to the distribution of load on the 
disk. This similarity is discussed further in part II. 

It should be noted that the distribution given for the 
uniformly loaded rotor is not the same as that given in 
reference 6, since in the present report the induced velocity 
is shown as being uniform across the lateral axis. A proof 
that the induced-velocity distribution must be uniform at this 
location is given in part II. In discussions between the 
senior author of reference 6 and the authors of the present 
report, it was suggested that the inaccuracy in the plane 
of the* disk resulted from a basic difficulty in the computa- 
tional procedure in the location near the edge of the disk. 
It is assumed that the values given in reference 6 are correct 
for locations outside the rotor tips. The calculations made 
for the present report were based on composite curves similar 
to the one shown in figure 13. 
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(b) x = 20.50° = tan -1 1/2. 
Figure 5. — Continued. 
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Figure 5. — Continued. 
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Figure 5. — Continued. 


INDUCED VELOCITIES NEAR A LIFTING ROTOR WITH NONUNIFORM DISK LOADING 



* 

N 


X/R 


REPORT 1319 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



430121 — 58 



(a) x = 63 . 43 ° = tan-i 2 . 

Figure 6.— Lines of constant value of induced velocity ratio v/vj in the longitudinal plane of symmetry of a rotor with a typical measured mean disk load.' 
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Figure 6. — Continued. 
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Figure 7. — Concluded. 
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Figure 8. — Induced velocity distribution along the X-axis of a rotor with a typical measured mean disk load. 
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±Y/R 

Figure 9. — Induced velocity distribution on the U-axis of a rotor with 
a triangular disk load. 




Figure 11. — Lines of constant value of induced velocity ratio v/vq in 
the longitudinal plane of symmetry of a rotor with a uniform disk 
load. x= 63.43° = tan -1 2. (Reproduced from ref. 6.) 



Figure 10. — Induced velocity distribution on the E-axis of a rotor 
with a typical measured mean disk loading. 


Figure 12. — Effect of disk load distribution on the distribution of 
induced velocity along X-axis. x = 63.43° = tan -1 2. 
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Fic.tre 13 . -Effect of disk load distribution on the distribution of 
induced velocity along the Y-axis. x = 90° = tan -1 


An interesting observation can be made from figure 13. 
Reference 1 lias shown that the uniformly loaded rotor values 
of induced velocity in the plane of symmetry (ref. 6) can be 
used as a guide to estimate the average induced velocity 
across a span of the order of the rotor diameter at any loca- 


tion. The downwash values for the lateral axis justify this 
result to the same order of accuracy as the measurements of 
reference 1, since the average induced velocity across the 
span for both nonuniform disk loadings is approximately the 
same as that for the uniformly loaded rotor. (That is, the 
average value of v/v 0 is about 1.0 in all three cases.) 
However, the differences between the results of reference 6 
and this report in the longitudinal plane of symmetry indi- 
cate the possibility of large errors if the uniformly loaded 
rotor flow field is used to estimate the effect of the rotor on 
items of appreciably different span than the rotor. 

CONCLUDING REMARKS FOR PART I 

A method of calculating the effect of nonuniform circu- 
larly symmetrical disk-load distributions on the normal 
component of the induced velocity of a lifting rotor has been 
presented in part 1. 

The induced velocity at the center of any rotor which has 
zero load at its center (including all practical rotors) must 
be zero. 

("harts have been presented for the normal component of 
induced velocity along the major axes and for the longi- 
tudinal plane of symmetry over a wide range of skew angles 
for rotors with two different, nonuniform axisymmetric disk 
loadings that are representative of the actual loading on the 
rotor disk. 

Rotor disk-load distribution has a large effect on the 
induced-velocity distribution and must be taken into account 
in estimating the effect of the rotor on most components of 
an aircraft. 


II. SYMMETRY RELATIONS AND THE RELATION BETWEEN RADIAL LOAD DISTRIBUTION 
AND THE RADIAL DOWNWASH DISTRIBUTION IN THE WAKE 

By S. Katzoff 


SYMMETRY RELATIONS CONCERNING THE INDUCED 
VELOCITIES IN THE PLANE OF THE ROTOR 

The first section of part II is concerned with the induced 
field of the skewed-cylinder vortex that is assumed in refer- 
ences 3 and 6 to represent the wake of a uniformly loaded 
disk. It will be shown, without reference to the detailed 
equations of the flow, that the induced-velocity field pro- 
duced in the plane of the disk by this idealized wake possesses 
certain symmetries. By providing relations between the 
induced velocities at pairs of symmetrically located points, 
such symmetries are useful for checking the accuracy of 
calculated values such as those of reference 6. They also 
provide directly the values of induced velocities along the 
lateral axis of the disk. 

As previously mentioned, the symbol notation followed in 
this part of the report differs from that usually used in rotary- 
wing work. The present notation was used in order to 
retain the usual complex-variable symbols in two-dimensional 
flow analyses and also to avoid a confusing multiplicity of 
subscripts. The reader should note in particular that the 
normal component of induced velocity' w is positive upward 
in the following material. 


POINTS ON THE KOTOR DISK 

In figure 14 (a) consider a pair of points on the rotor disk, 
such as P and Q , that have the same lateral location and are 
symmetrically located with respect to the lateral axis. That 
is, if P is at point (#, y, 0), then Q is at point (— x, y, 0). 
Let the three components of the induced velocity at P be 
(u\y V \ , w x ) and at Q be ( u 2 , v 2 , w 2 ), where the third component 
is generally the only component of interest. 

Now consider the entire diagram rotated 180° about the 
F-axis so that it appears as in figure 14 (b). Points P and 
Q thereby move to points P' and Q' , where the induced 
velocities are, respectively, (— U\, v u —w x ) and (— u 2 , v 2 , 
—w 2 ). 

If figure 14 (b), with the vorticity reversed, is fitted to 
figure 14 (a), so that P' falls on Q, and Q' falls on P, the two 
semi-infinite vortex cylinders form a single continuous 
infinite vortex cylinder (fig. 14 (c)). The induced velocity 
at the left-hand point (P 7 and Q) is ( u x -\-u 2 , v 2 —v u Wi-\-w 2 ), 
and at the right-hand point (Q' and P) it is (iq-+-w 2 , Vi~v 2 , 
Wi+w 2 ). The two velocities are thus equal except for the 
lateral components, which are equal and opposite. 

In appendix A of reference 3, however, it was shown that 
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(a) Rotor and wake. (b) Wake rotated 180° about y-axis. 

(c) Vortex cylinders shown in figures 14 (a) and (b) fitted together to 
form an infinite cylinder. 

Figure 14. — Symmetrically located points in the rotor disk. 

the induced velocity within this infinite skewed helix is 
uniform and has no lateral component. The first result of 
the present discussion, then, is that v 2 —Vi=Vi~v 2 =0 or 
Vi=V 2 . 

The second result is that the sums Ui~\~u 2 and Wi~\~w 2 are 
uniform. That is, the sums of the two longitudinal com- 
ponents and of the two vertical components are the same 
for all pairs of symmetrically located points on the disk (as 
P and Q), and are equal to the longitudinal and vertical 
components of the induced velocity within the infinitely 
long helix. In the nomenclature of reference 3, 

Ui~\-U 2 =V x ' 

Wi +w 2 =v/ 

Furthermore, if points P and Q are on the lateral (IQ axis, 
so that they coincide (u { = u 2 and Wi=w 2 ), 



That is, the longitudinal and vertical components of the 
induced velocity are constant all along the lateral diameter 
of the disk and are equal to the values at the center of the 
disk. 

POINTS OUTSIDE THE ROTOR DISK 

If points P and Q are beyond the edge of the rotor disk 
(fig. 15), the discussion proceeds as before except that it no 
longer follows that u x -\-u 2 and w x -\-w 2 are constant or that 
Vi—v 2 is zero. The reason, of course, is that after super- 
position of the two semi-infinite skewed vortex cylinders 


z 



Figure 15. — Symmetrically located points in the plane of, but outside, 

the rotor disk. 

the points lie outside the infinite elliptic cylinder instead of 
inside of it. The combined flow field outside the cylinder 
is not uniform but is the same as if the cylinder had a solid 
boundary. That is, the combined induced flow is two- 
dimensional in planes normal to the cylinder axis and can 
be computed as that induced by an elliptic cylinder in a 
steady flow of velocity V cos x normal to its axis. The 
streamlines of the induced flow in the normal plane are 
sketched in figure 16. The components u' and if of the 
induced flow in this plane are related to Ui+ik, Vi—v 2 , and 
Wi~\~w 2 as follows: 

Ui~\~U 2 = v' COS X 
w 1 J rW 2 =v' sin X 

V\ — v 2 =u' at the right-hand point (P and Q') 

For points along the lateral axis, 

Ui-\-u 2 =2ui=v' cos x 
w l J r w 2 =2w l =v / sin X 

Vi—v 2 =0 

Thus, for points along the lateral axis, u x and w x can be 
determined directly 7 from the known two-dimensional flow 
about the ellipse. 
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Figure 16. — Two-dimensional crossflow around an ellipse. 
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In order to facilitate application of these results, a few 
remarks concerning the calculation of v' are contained in the 
next section. 


FLOW ABOUT AN ELLIPTIC CYLINDER 

Consider (fig. 16) the flow induced by an elliptic cylinder 
of semimajor axis a and semiminor axis b in a cross flow of 
velocity V cos x- Let the plane of the figure be the z' plane, 
where z'=x'-\-iy ' (the x ' axis is parallel to the y- axis, or the 
lateral axis of the disk). 

The flow is probably best treated in elliptic coordinates, 
related to the rectangular coordinates by 


where 


z'=c cosh f 
c =JcL^b 2 


The complex flow function w for the field induced by the 
ellipse in the steady cross flow of velocity V cos x is 

w=(V cos x) 

(The equation may be obtained, for example, from the 
equation on page 256 in ref. 8 by setting U equal to zero, 
replacing V by V cos x , and omitting the free-stream flow 
function.) Then the complex velocity is 


dw 

dz' 


dw/d£ 
dz '/d£ 


=—(V cos X) ia 


v 


a+b e f 
a—bc sinh 


f 


—{V cos x)^ 


l a+b J2_ 
a— b e 2 *— 1 


Setting f =£ +it] and multiplying numerator and denominator 
by the conjugate of the denominator (e 2 ^~ ir,) — \) leads to 
the following expression for v', the negative of the imaginary 

* r dw 
part of 

dz 


y'=2(F cos X) 


u 


a-\-b e 2 * cos 2tj— 1 
a—be 4 i—2e 2 i cos 2r/+l 


or 


v'=(V cos x) 


a(a-\-b) cos 2r) — e 2 * 
c 2 cosh 2£— cos 2rj 


Since x'+iy'=c cosh (£+irj), the values of £ and rj correspond- 
ing to any point (x ' , y') in the field of the ellipse are obtain- 
able directly from tables of complex hyperbolic functions 
(ref. 9). (In using these tables, note that rj is given in units 

b . c \ 

of 7r/2. It will be observed that cos and sin X=— ) 

a a / 

The points P, P', Q, and Q' will be noted to lie in a plane 
which intersects the plane of the ellipse at the angle x- 
However, because of the two-dimensional character of the 
flow, the desired values of the induced velocities at a point 


( x , y) in the plane of the rotor may be found at the point 
(x' = y, y'=x cos x) in the plane of the ellipse. 

RELATIONS BETWEEN THE RADIAL VARIATION OF DISK 
LOADING AND THE RADIAL VARIATION OF DOWNWASH 
VELOCITY WITHIN THE DOWNSTREAM WAKE 

In figure 13 the computed lateral variation of the down- 
wash angle across the lateral axis was shown to resemble 
the assumed radial distribution of disk loading. An effort 
was therefore made to determine whether the resemblance 
was mainly fortuitous or whether a general theoretical basis 
existed for it. 

It was found possible to show that, for the far wake, the 
downwash velocity along any radius is proportional to r n if 
the disk loading is also proportional to r n . The propor- 
tionality factor, however, depends on n (in addition to the 
azimuth). Accordingly, if the radial disk loading is repre- 
sented by a power series, the downwash velocity along any 
radius is also represented by a power series where, however, 
the respective coefficients of the two series are not propor- 
tional. Thus the radial loading is not, in general, exactly 
proportional to the radial downwash-velocity distribution, 
although the two may be nearly proportional if the radial 
loading is approximately proportional to, say, r or r 2 . 

In the following development the simplest case, x=0, 
will be discussed first. The general case, 0 o< Cx^90°, will 
then be discussed. Finally, the particular case x=9A°, 
which is related to linearized lifting-surface theory, will be 
discussed. 

THE CASE x = 0° 

The problem is simple for the circular wake, x=0°- 
Consider first the uniformly loaded disk, which produces a 
single helical vortex lying along a circular cylinder. In the 
far wake the downwash within the cylinder is uniform, and 
there is no downwash at all outside the cylinder. For an 
infinitesimal annulus (dr) of disk loading, the wake is then 
two concentric circular cylinders having equal and opposite 
vorticity and separated by dr. The downwash between the 
two cylinders is uniform and proportional to the annulus 
loading. The downwash everywhere else is zero. If the 
disk loading is made up of a continuous distribution of such 
annuli, the downwash directly downstream of each annulus is 
thus determined only by that annulus loading and is not 
affected by the loading of any other annulus. For this case, 
then, the radial variation of downwash velocity is propor- 
tional to the radial variation of disk loading. Restricting 
the theorem to loadings that are proportional to r n is unneces- 
sary for this case, which will be recognized as merely the well- 
known propeller blade-element theory applied to the far 
wake instead of to the disk itself where the induced velocities 
are only half as much. 

THE CASE 0°<x<90° 

As previously mentioned, the uniformly loaded disk pro- 
duces an ultimate wake in the form of an inclined elliptic 
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cylinder, within which the induced flow is uniform and outside 
of which the induced flow is that due to an elliptic cylinder in 
the cross-flow component of the free-stream flow. If the 
loading is circularly symmetrical but not uniform, the wake 
consists of a corresponding distribution of concentric similar 
elliptic cylinders, and the present problem is concerned with 
the superposition of their induced flow fields. 

If each cylinder induced only the uniform internal velocity 
field and induced no external field, the argument would 
proceed just as for the case x=0°, and it would follow im- 
mediately that the radial distribution of induced velocity is 
proportional to the radial loading distribution. In order to 
prove the theorem, then, it is only necessary to show that the 
total of the external fields of the distribution of cylinders will, 
along any radius, also have velocities proportional to r n , if 
the loading is proportional to r n . 

Consider the flow’ induced by an ellipse of given shape and 
unit size (say unit semimajor axis) moving at velocity V cos x. 
As previously indicated, the induced velocity component 
parallel to the minor axis (parallel to the direction of motion 
of the ellipse) is the component of interest. Along any di- 
rection 0 from the origin, this induced velocity is a function 
of the radial distance p. 

If, while the shape remains constant, the linear dimensions 
of the ellipse are changed by the factor k, the velocity 
originally existing at radial distance p will be found at 
radial distance kp along the given direction 0. For the 
general ellipse of size k } then, the velocity component v f can 
be expressed as 


v'=Mk/ P ) 

The velocity induced by an increment in size, or annulus, 
corresponding to dk is 


to p*, which was to be proved. The proportionality factor 
is a function of both 0 and n. It will be observed t hat this 
part of the proof can apply to any shape, since the fact that 
the wake is elliptical was not required or used in the proof. 

Performing the indicated integration in order to get the 
proportionality factor F(6,n) is obviously a somewhat 
awkward task. Possibly the most convenient method is 
first to integrate by parts in order to return /' to / (that is, 
v ' , for which an expression was previously presented), and 
then determine the resulting integral numerically. For 
either the lateral or longitudinal axis, however, the calcula- 
tion may be performed without excessive difficulty. The 
method will be here indicated. 

For the lateral axis ( 77 = 0 ), the previously derived expres- 
sion for v' reduces to 


v'=(V cos x) 


a /a+ 
c \ a— 


+ b x'-m'x'I-c 


V* 


so that the problem is reduced to evaluating 


r° x \ n t> ( x'-^' 2 -c 2 \ 

Jc-o M / 


dc 


Substituting c=x' sin a changes the integral to a form that 
can be integrated by application of items 274 and 263 of 
reference 10. 

For the longitudinal axis the expression for v' reduces 
to 




d ^ d k=f(d,k/p) <l p k 

where /' denotes the derivative with respect to k/p. If each 
incremental velocity is weighted by the factor k n , the total 
velocity is 

r k”f(e,k/p)-= P ” r (* /pyf'(o,k/p)~ 

Jk = 0 P Jk= 0 ‘ p 

where the upper limit K is that value of k for which the ellipse 
just touches the specified point (p,0). (As previously noted, 
the present discussion concerns only the contributions of 
those annuli for which the point is located externally.) 
Changing the variable of integration from k to Ar/p changes 
the form of the preceding expression to 

J 'K/p 

(k/pTf ( e,k/ p )d(k/p ) = p”F(e,n) 

k/p = 0 

since the upper limit K/p is now a function only of 0. Thus, 
the induced velocity along any radial line is proportional 


so t hat the problem is reduced to evaluating 


p v c „ a fffiY) 

J c=0 <*\ Vc 2 +j /' 2 / 


dc 


Substituting c=y f tan a changes the integral to a form that 
can be evaluated by application of item 274 of reference 10. 


THE CASE x = 90° 

The limiting case of the flat wake is frequently assumed 
for convenience in analyses pertaining to high-speed forward 
flight. Since the rotor and the wake are in the same plane 
the case is analogous to that of classical linearized wing 
theory, with the wing in this instance having a circular plan 
form and a circularly symmetrical loading. 

Consider first a uniformly loaded rotor. The “span load” 
distribution of the disk, considered as a wing, is then ellip- 
tical, since it is merely proportional, at each spanwise station, 
to the “chord.” As is well known, the down wash in the 
far wake behind an elliptically loaded wing is uniform. That 
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is, disregarding the rolling-up phenomenon, the flow in the 
vertical plane normal to the far wake is the two-dimensional 
flow about a horizontal straight line of span equal to the 
rotor diameter, moving downward with a velocity equal to 


2 r L v 


7r (Aspect ratio) 


o p T T 

= - 7a l (see fig. 17). The complex ve- 

7r(4/7r) ° 


C T \ 


locitv for this two-dimensional flow around the 
given by 





+ i 


far 


wake 


is 


where 

w complex flow function 

r semispan of wing (or radius of the rotor) 

z'=x'-\-iy' where, in the present case (x = 90°), x' is the 
same as lateral dimension ?/ and y ' is the same 
as vertical dimension z 

If the radius of the uniformly loaded circular wing is 
increased by dr, the complex veocity in the flow field changes 
by 


d (dv)\j C L V\ -j irz ' 
Tr\dz ') ar ~ 2 L(^ ,2 -^ 2 ) 3/2 J 


dr 


This expression thus gives the complex velocity in the far 
wake contributed by a wing in the form of a circular annulus 
of radius r and width dr, and having a lift coefficient C L . 

Now consider a circular wing on which the loading is 
proportional to the ?ith power of the radius. It is made up 
of annuli of variable radius r and lift coefficient given by, 
say, Ar n , where A is a constant. The total complex velocity 
in the far wake is found by substituting Ar n for C L in the 
preceding expression and integrating with respect to r 
between 0 and the outermost annulus, r—R. The complex 
velocity is thus 



where the upper limit R is greater than |z'|. 



Figure 17. — Two-dimensional crossflow around a straight line. 


In order to avoid a singular point in the integration, the 
point z r will be assumed to lie slightly above the real axis 
and will be allowed to approach the real axis after the 
integration is performed. Substituting r=z' sin d (where d, 
in general, is complex) transforms the integral to 

— iAz’ n V f sin sm n+1 6dd 


\z' n V 2 sin w+1 0 
2 J* =0 cos 2 1 


which can be evaluated by first applying the second form 
of item 274 of reference 10 and then reducing y'sin” +1 0 dd by 
successive applications of item 263. Except for the last 
term, the integral is 

— i4z' n F^sin n+1 d tan0+sin n 0 cos 0-f^ - sin”" 2 0 cos 0+ 

r ~ iV ~ sinn ~ 4 6 c ° s g + • • 1 

(n— l)(n— 3) J 

Substituting the limits and lettering z' become a pure real, 
less than R, results in a pure real; that is, all of these terms 
contribute nothing to the vertical-velocity component. The 
final term of the integral involves 


-i- , -i" 

Csin 2 /*sin 2 

sin 6 dd or 

Jo Jo 


dd 


depending on whether n is even or odd, respectively. For 
n even, 


R. I 


. -i4 

sin 2 


Csin 2' 

3 . sin d dd=R. P.(— cos d) 

=R.P. [-V>-S+l] 

which approaches 1 as z’ approaches a pure real, x’ < B, 
since the first term in the bracket approaches a pure imagi- 
nary. For n odd, 


R. P 


j: 


. -i* 
s\n 2 ' 


de= R. P.( 0 ) 


which may be evaluated as follows: 

Let 

6=p+iq 

then 

sin 0=sin (p-f ig)=sin j> cosh g+i cos p sinli <}= . , 

n 

Since approaches a pure real as the imaginary 

term must be zero; that, is, cos p= 0, or p = 7 r/ 2 . Accordingly, 


R. P. (0) 


. -i* 
sin 2 ' 


>7t/2 


as z' approaches a pure real, x'<^R. The final result for 
the downwash velocity at the real point z' —x' is, for n odd, 

Ax' n V T n(n—2)(n— 4) . . . (1) "V 
2 L (n— 1)(^— 3) . . .(2) J2 
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and, for n even, 



n(n-2)(n-4) . . . (2)"] 

.(ft— 1) (ft— 3) . . . (1) J 


For n= 0, the result is -^-5 for n= 1, it is ^ ' ^ 5 for n= 2, 
it is VAx /2 ; for ?i— 3, it is ^ 7rF^4a:' 3 . 

o 


CONCLUSIONS FOR PART II 

A study of the symmetry relations concerning the induced 
velocities in the plane of a uniformly loaded rotor indicates 
that: 

1. The sum of the induced velocities at points (x,y) and 


(—x,y) in the rotor disk is constant and equal to twice the 
induced velocity at the center. 

2. The induced velocity is constant all along the lateral 
diameter of the disk. 

3. The sum of the induced velocities at points (x } y) and 
(—x,y) in the plane of, but outside, the rotor disk equals the 
induced velocity at the corresponding point near the far 
wake. This sum can be determined from the two-dimensional 
flow about the ellipse that represents the wake cross section. 

A study of the relation between the radial load distribution 
and the radial downwash distribution in the wake of a non- 
uniformly loaded rotor shows that, if the disk § loading varies 
as the nth power of the radius, the induced velocity in the 
far wake also varies as the nth power of the radius. 


III. ANALYSIS AND COMPARISON WITH THEORY OF FLOW-FIELD MEASUREMENTS NEAR A 
LIFTING ROTOR IN THE LANGLEY FULL-SCALE TUNNEL 

By Harry H. Heyson 


APPARATUS AND TESTS 

The surveys wtTxT^ohdticted in the Langley full-scale 
tunnel, which is described in reference 1 1. The experimental 
apparatus used in these tests is shown in figure 18 (a). 

The rotor was of the teetering or seesaw type and had 
un tapered, untwisted blades (fig. 18 (b)) with NACA 0012 
airfoil section. The rotor radius was 7.5 feet and its solidity 
was 0.0543. The rotor tip speed was 500 feet per second for 
all tests except that at n= 0.232. For this test, the rotor 
was operated at 450 feet per second. 

Two different survey rakes were used in the tests (fig. 
18 (c)). A calibrated five-tube survey rake (fig. 18 (d)) was 
used to measure the stream pitch and yaw angles and 
dynamic pressures for most of the tests. This rake is fully 



described in reference 12. In order to obtain data closer to 
the plane of the rotor, a calibrated pitch head (fig. 18 (e)) 
was used to survey points in the longitudinal plane of sym- 
metry tor one flight condition. 

The procedure used in the tests was as follows: The rotor 
was firsbiset approximately at a predetermined flight condi- 
Then, the area above, below, and behind the rotor 


tlOIL 


was surveyed with the rakes. At intervals during the 
surveying, readings were taken of rotor thrust and drag, 
blade flapping and feathering motions, and tunnel dynamic 
pressure. 


, 


CORRECTIONS 


A floor was installed in the tunnel test section for all the 
tests, and for t his configuration the Langley full-scale tunnel 
has no appreciable jet-boundary correction (ref. 13). r 
In order to correct the measured dynamic pressures and 
stream angles- for the flow about the supporting mechanism 
and for tunnel stream angle, data were taken twice at each 
point in space — once with and once without the rotor blades 
installed. The difference is considered to be characteristic 



(a) Test setup looking into entrance cone of tunnel. 
Fic.ure 18. — -Equipment used in surveys. 


(b) Rotor blade. 
Figcre 18. — Continued. 
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(c) Rakes as mounted on survey-carriage ram. 
Figure 18. — Continued. 



(d) End of one tube of five-tube rake. 
Figure 18. — Continued. 


of the isolated rotor. It should be noted, however, that this 
procedure is not adequate in regions immediately behind 
the rotor. Measurements made in these regions have little 
accuracy because of the effect of the rotor on the flow about 
the tower. The angle of attack of the rotor was corrected 
for tunnel stream angle. 


PRECISION 

The accuracy' of the two rakes is essentially the same. 
Both measure dynamic pressure within 2 percent and stream 
angles within y 2 °. These figures take into account the un- 
steadiness of the flow and the high vibration level of the 
rake mount. At the lowest tip-speed ratio, additional in- 
accuracy is caused by the difficulty in reading the small 
changes of liquid height in the manometer tubes. Overall 



(e) Small pitch head. 
Figure 18. —Concluded. 


accuracy in terms of the average value of the normal com- 
ponent of induced velocity will be a function of the rotor lift 
coefficient. For ease of reference, the following table gives 
the estimated overall accuracy for each of the flight condi- 
tions covered in this investigation: 


Tip-speed ratio, m 

Accuracy, percent v 0 

0. 095 

±15 

. 139 

±10 

. 140 

±10 

. 232 

±25 


The rakes used were incapable of measuring instantaneous 
variations of induced velocity such as might be desired for a 
refined blade-vibration analysis. Therefore, all data are 
time-averaged values. The experimental results of refer- 
ence 2 indicate that instantaneous values may differ greatly 
from the mean values measured in these tests. 

During the tests, the rotor tip speed was held constant to 
within 1 percent. The thrust coefficient was measured within 
1 percent. The skew angle, which was computed by using 
the measured values of C T and the measured blade motions 
in equation (29) of reference 6, should be accurate to within 1°. 

The position of the rake with respect to the hub was known 
within 0.3 percent R. Because the rotor was of the teetering 
type, the only coning was due to the deflection of the blades. 
The tip deflection was estimated visually to be between 2 
and 3 percent of the radius. The data are presented with 
reference to a plane parallel to the tip-path plane and passing 
through the teetering pin. 
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RESULTS AND DISCUSSION 
BASIC DATA 

Because of the present lack of information, both theoretical 
and experimental, on the induced velocities in many areas 
near the rotor, substantially all the basic data obtained in 
the investigation are presented in the appendix. In each case 
the rotor flight condition is identified on the figure by the 
skew angle X and the tip-speed ratio /z. Table I identifies the 
various flight conditions more completely. 

The basic-data curves are not discussed as such in this 
report since the more significant phenomena have already 
been pointed out in reference 1. 

DISCUSSION OF APPLICABLE THEORY 

In reference 6, Castles and De Leeuw present calculated 
values of the normal component of induced velocity in the 
longitudinal plane of symmetry and along the principal 
(A", F,Z) axes of a rotor with a uniform disk loading. Part I 
uses the flow field of reference 6 to calculate the flow field 
for rotors having non uniform but circularly symmetrical disk 
loadings (the specific cases calculated being rotors with a 
triangular disk loading and with a typical measured mean 
disk loading). 

Additional theoretical calculations by Mangier and Squire 
(ref. 5) give induced- velocity information in the entire rotor 
disk and also for a transverse plane in the far wake. A later 
section discusses a method of transforming these results to 
correspond with the results of reference 6 and part I except 
for the specific disk loadings considered. 

Drees (ref. 4) discusses some of the effects of a finite tip- 
speed ratio and presents calculations of some of these effects 
at selected points on the rotor disk. It should be noted that 
the calculated curves presented with the data are not modi- 
fied for the asymmetries discussed in this reference. 

COMPARISON OF MEASURED DATA WITH CALCULATED 
VALUES OF REFERENCE 6 AND PART I 

Longitudinal plane of symmetry. — Figures 19 to 22 present 
the measured values of induced velocity in the longitudinal 
plane of symmetry as compared with the calculated values 
from reference 6 and part I. 

For the four flight conditions, the measured data near the 
forward portions of the rotor fall close to the flow calculated 
for the rotors with zero load in the center. The existence of 
essentially zero induced velocity at the center of the rotor 
and of regions of upwash behind the center are clearly shown. 
Further rearward in the flow, just behind the trailing edge of 
the rotor {X/R = 1.07), there is no longer any pronounced 
correlation between the measured and the calculated flow, 


regardless of the lift distribution assumed to exist on 
the rotor disk. At greater distances behind the rotor 
{XI R =2. 01 and 3.14), there seems to be no correlation. 
The physical reasons for this lack of correlation are discussed 
subsequently. 

Longitudinal axis. — If curves are faired through the data 
in figures 19 to 22 and the values at Z/R = 0 are taken from 
these curves, it is possible to obtain an idea of the variation 
of induced velocity along the longitudinal axis of the rotor. 
The results, plotted in figure 23, show clearly the good corre- 
lation of the nonuniformly loaded rotor calculations for the 
flow over the forward three-quarters of the rotor and the 
deviation behind that point. 

Figure 24 compares the measured points of figures 23 (c) 
and 23 (d). These two flight conditions have very small 
differences in skew angles. There is, however, an appreciable 
difference in tip-speed ratio and lift coefficient. Actually, 
the calculated induced-velocity ratios for these two skew 
angles are virtually identical. It may be seen that, experi- 
mentally, there is very little difference between the induced 
velocities along the longitudinal axis for the two tip-speed 
ratios. This fact is as assumed in references 3 and 6. 

Lateral axis. The basic data may be cross-plotted in the 
same manner to obtain the induced-velocity distribution 
along the lateral (F) axis. The results of these cross plots 
are shown in figure 25 and are compared with the calculated 
distributions of reference 6 and part I. The calculated flow 
is symmetrical since no account is taken of dissymmetries 
due to the finite tip-speed ratio. 

A word of caution is necessary at this point. Because of 
the necessity of maintaining a clearance between the survey 
rake and the rotor, it was not possible to obtain data much 
closer than 20 percent R above the rotor. (See figs. 19, 2 1 
and 22.) The theory predicts rather large gradients of 
induced velocity at locations closer to the rotor than this 
point. In the longitudinal plane of symmetry it is possible 
to cross-plot the data with some confidence by comparing 
the trends shown well above the rotor by the theory and by 
the measured points. At other locations, without the theory 
as a guide, the accuracy of cross plots is problematical. 
Therefore, the following discussion of the induced- velocity 
distribution on the F-axis should be regarded as being on a 
less firm basis than the discussion in the preceding sections. 

Figure 25 indicates that the rotor does have zero induced 
velocity at its center. This is in agreement with the assump- 
tion that the rotor carries zero load at its center and further 
indicates that this fact must be used if the flow is to be 
calculated with any degree of accuracy. 


TABLE I.— FLIGHT CONDITIONS 


Tip-speed 
ratio, n 

Rotor wake 
skew angle, 
X, deg 

C L 

Ct 

Mi, 

ft/sec 

Disk loading, 
lb/sq ft 

Equivalent flat- 
plate area, /, 
sq ft 

Angle of 
attack, 
ot, deg 

0. 095 

75. 0 

0. 720 

0. 00320 

500 

1. 91 

2. 77 

-9. 2 

. 13 9 

85. 8 

380 

00373 

500 

*2. 18 

2. 34 

- 1 . 1 

. 140 

82. 3 

373 

«. 00371 

500 

a 2. 21 

2. 79 

-5. 3 

. 232 

83. 9 

. 122 

. 00321 

450 

1. 55 

1. 35 

-9. 5 


® The seeming inconsistency between the values of rotor lift coefficient and thrust coefficient on one hand and the values of average disk load on the other wa 
caused by a slight difference in the density of air at the different times during which the tests were run. 
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(a) X/R=- 1.0. (b) X/R=- 0.5. (c) X/R= 0. (d) X/R = 0.5. 

(e) X/R= 1.07. (f) X/R= 2.07. (g) XIR = 3.U. 

Figure 19. — Comparison of measured and theoretical values of induced-velocity ratio v/vo in longitudinal plane of symmetry. x = ~5.()°; ai=0.095. 
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Z/R 



(a) X/R=- 0.95. (b) X/R = 0.8. (c) X/R=- 0.6. (d) X/R 0.4. 

(e) X/R= —0.2. (f) X/R- 0. (g) X/R -0.2. 

(h) X/R-0A. (i) X/R = 0.0. 

(j) X/R=0.8. (k) X/R — 0.05. 

Figure 20. Comparison of measured and theoretical values of induced- velocity ratio v/tio in longitudinal plane of symmetry, x = 85.8°; ^ 0.139. 


•Ui(M2l -58- 
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(a) X/R = — 1.0. (b) X/tf=- 0.5. (c) X/R = 0. (d) X//? = 0.5. 

(e) X/«=1.07. (f) A'//? = 2.07. (g) X/fl=3.14. 

Figure 21. — Comparison of measured and theoretical values of induced-velocity ratio v/v a in longitudinal plane of symmetry. x = S2.3°; ^ = 0.140. 
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(a) X/R=- 1.0. (b) X/R=- 0.5. (c) X/R = 0. (d) X/R = 0.5. 

(e) X/R=l.07. ({) X/R = 2.07. (g) X/R = 3.14. 

Figure 22. Comparison of measured and theoretical values of induced-velocity ratio v/v 0 in longitudinal plane of symmetry. x=83.i)°; ^ = 0.232. 
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(a) x = "5.0° ; M = 0.095. 

(b) *=85.8°; M =0. 139. 

(c) x=82.3°; m = (U4(). 

(d) x =83.9°; m=0.232. 

Figure 23.- -Comparison of measured and theoretical values of induced-velocity ratio v/v 0 along longitudinal center line. 
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Figure 24. — Effect of tip-speed ratio on distribution of induced- 
velocity ratio v/v 0 along longitudinal center line. 



(a) x= 75.0° ; /x= 0.095. 

(b) x = 82.3° ; 0.1 40. 

(c) x = 83.9°; /z = 0.232. 

Figure 25.— Comparison of measured and theoretical values of induced- 
velocity ratio v/vq along lateral center line. 

Figure 25 also shows that the induced- velocity peak on 
the retreating side is higher and further outboard than the 
peak on the advancing side. This result is logical when the 
difference in loading on the two sides of the rotor are com- 
pared. The resultant velocities at the inner portion of the 
blades on the retreating side are very low or even negative, 
whereas the corresponding velocities on the advancing side 


are much higher. Therefore, the loading on the inner por- 
tion of the retreating side of the disk must be less than that 
on the advancing side. The outer portion of the retreating 
side must then carry a greater loading in order that the blade 
thrust moment remain essentially constant around the disk. 
Thus, the loading on the retreating side will peak farther out- 
board and at a higher value than the loading on the advanc- 
ing side of the disk. This is the same trend shown by the 
measured induced velocities and is indicative of the very 
close relation between the induced velocities and the loading 
along this axis. 


COMPARISON OF MEASURED DATA WITH CALCULATED VALUES 
OF REFERENCE 5 

Modification of calculations for finite lift coefficient. — 

Mangier and Squire, in reference 5, compute the induced 
velocities in the rotor disk and in a transverse plane in the 
far wake. Their calculated results are limited to vanishingly 
small lift coefficients in that the trailing-vortex system leaves 
the rotor in the same direction and with the same velocity 
as the free stream. 

The assumption of a finite lift coefficient alters this picture 
in two ways: First, the trailing-vortex system must now 

cross the free-stream flow and, second, there will be a change 
in the vortex density. 

The results of reference 5 are easily modified to account 
for these changes and, when so modified, will correspond 
with the results of reference 6 and part I except for the 
specific disk loadings considered. The first modification is 
merely a rotation of the field to account for the inclination 
of the wake. The incidence i is then the complement of 
the skew angle X. The change in vortex density, although 
affecting the absolute induced velocity v , w ill have no effect 


upon the induced-velocity ratio 


vo 4 (ucj‘ 


Therefore, 


the new value of v may be found by using the value of v/v 0 
from reference 5 and by calculating v 0 from an equation 
which considers the finite lift coefficient (such as eq. (32) 
of ref. 6). 

For the p impose of comparison with the measured data 
of this report, all the induced velocities of reference 5 have 
been modified to correspond with reference 6 in regard to 
both skew r angle and vortex density. 

Plane of the rotor. — Values of induced-velocity ratio in the 
plane of the rotor obtained by cross-plotting the measured 
data are compared with the corrected calculated values from 
reference 5 in figures 26 to 32. These figures are discussed 
only in general, since they indicate the same results as the 
previous comparison w T ith reference 6 and part I. This re- 
sult might be expected since the calculated results should 
be comparable if the same load distributions are assumed. 
The resemblance between the typical loading of part I and 
pressure distribution III of reference 5 is evident. 

In general, the same extent of agreement is noticeable 
(sec especially fig. 26 (b)); that is, the calculated flow^ field 
is close to the measured field only over the forward three- 
quarters of the disk. 

Each figure (figs. 26 to 32) shows that the induced- 
velocity field for pressure distribution III is closer to the 
measured field on both sides of the disk than that for pres- 
sure distribution I. Since pressure distribution III repre- 
sents a disk loading with zero load in the center and the 
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(a) x=75.0°; m = 0.095. (b) x=85.5°; m = 0.139. 

(c) x=82.3°; m = 0.140. (d) x =83.9°; M = 0.232. 

Figure 26— Comparison of measured values of induced-velocity ratio v/v 0 on longitudinal center line with theory of reference 5 modified. 


other distribution does not, this result is as expected. 

This result does not agree with that of reference 14, which 
indicates that pressure distribution I should be used for 
calculating the induced velocities over the advancing side 
and pressure distribution III, over the retreating side of 
the rotor. This result (ref. 14) arose from measurements 
made at X/R=1.5 and, as such, is subject to error because 
of the discrepancy between theory and the measured data 
in this region. The present result is considered to be more 
accurate. 

The values at XIR= 0.5 (fig. 32) seem to indicate that, 
at this location, the error is generally somewhat larger near 
the outer edges of the disk. This result might be expected 


since the roll-up of the t railing-vortex system would be 
expected to affect the outer edges of the disk sooner than the 
center. 

Far wake. — Reference 5 also computes the induced-velocity 
distribution in the far wake. The' computed contours of 
induced-velocity ratio for both pressure distributions in the 
far wake are compared in figures 33 to 38 with measured 
contours of induced velocity at XjR=2.01 and X/R= 3.14 
for three different flight conditions. (Notice that the rota- 
tion of the calculated flow field required to obtain correspond- 
ence with reference 6 results in a vertical displacement of 
the computed induced velocities.) 
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X/R 


(a) x = 75.0°; *<=0.095. 

(b) x=82.3°; *<=0.140. 

(c) x=83.9°; **=0.232. 

Figure 27. — Comparison of measured values of induced-velocity ratio 
v/vo at Y/R = ±0.25 with theory of reference 5 modified. 



(a) x=75.0°; *<=0.095. 

(b) x=82.3°; **=0.140. (c) x=83.9°; **=0.232. 

Figure 28. — Comparison of measured values of induced-velocity ratio 
tt/tto at Y/R— ±0.50 with theory of reference 5 modified. 
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(a) x=75.0°; *. = 0.095. 

(b) x = 82.3° ; m = 0.140. (c) x = 83.9°; M = 0.232. 

Figure 29. Comparison of measured values of induced-velocity ratio v/v« at YIR=± 0.75 with theory of reference 5 modified. 
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(a) x=75.0°; M = 0.095. 

(b) x— 82.3°; ju=0.140. 

(c) x= 83.9° ; M = 0.232. 

Figure 30. — Comparison of measured values of induced-velocity ralio 
vlvi i at X/R-— 0.5 with theory of reference 5 modified. 



(a) x=75.0°; M =0.095. 

(b) x=82.3°; m =0.140. 

(c) x = 83.9°; /»=0.232. 

Figure 31. — Comparison of measured values of induced velocity ratio 
v/vo on lateral center line with theory of reference 5 modified. 


400421 — 08 G 
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Figure 32. — Comparison of measured values of induced-velocity ratio v/vt, at X/R=0.5 with theory of reference 5 modified. 
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Z/R 


- 1.0 



-1.5 


(a) Measured data. (b) Theoretical, pressure distribution I. (c) Theoretical, pressure distribution III. 

Figure 33.— Comparison of measured values of induced-velocity ratio v/v 0 at X/R= 2.07 with theory of reference 5 modified. x =75.0° ; „ 


= 0.095. 



40 


REPORT 1319 NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Advancing side Y/R Retreating side 



1.0 .5 0 -.5 -1.0 -1.5 


Y/R 



(a) Measured data. 

Figure 34. — Comparison of measured 


(b) Theoretical, pressure distribution I. (c) Theoretical, pressure distribution III. 

values of induced-velocity ratio c/r 0 at X//f = 3.14 with theory of reference 5 modified. x=75.0°; 


= 0.09 
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Z/R 




Y / R 


Z/R 


- 1.0 



- 1.5 


(a) Measured data. (b) Theoretical, pressure distribution I. (c) Theoretical, pressure distribution III. 

Figure 35— Comparison of measured values of induced-velocity ratio v/v 0 at X/R = 2.07 with theory of reference 5 modified, x=82.3°; , 


430421 — 58 7 


= 0.140. 
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Z/R 


-.5 


- 1.0 



-1.5 


(a) Measured data (b) Theoretical, pressure distribution I. (c) Theoretical, pressure distribution III. 

Figure 36 — Comparison of measured values of induced-velocity ratio v/v 0 at X/R = 3.14 with theory of reference 5 modified. x=82.3°; n 


= 0.140. 
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Advancing side 


Y/R 


Retreating side 



1.0 1 1 ^ * 

1.5 1.0 .5 0 -.5 -1.0 -1.5 


Y/R 


Z/R 



1.5 


I .0 


.5 


0 

Y/R 


-.5 


- 1.0 


-1.5 


(a) Measured data (b) Theoretical., pressure distribution I. (c) Theoretical, pressure distribution III. 

i (jure 37. — Comparison of measured values of induced- velocity ratio v/v 0 at X/R = 2.07 with theory of reference 5 modified. x=83.0°; n 


= 0.232 
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1.5 


I .0 


.5 


0 

Y/R 


- .5 


- 1.0 


-1.5 


(a) Measured data. (b) Theoretical, pressure distribution I. (c) Theoretical, pressure distribution III. 

Figure 38. — Comparison of measured values of induced-velocity ratio v/v 0 at A//? = 3. 14 with theory of reference 5 modified, x— 83.9°; / 


= 0.232. 
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In general, the measured induced velocities are of the 
same order of magnitude as the calculated values. This 
result would be expected, regardless of the flow pattern 
assumed for the calculations, since a given amount of 
downward momentum must be imparted to the air to 
obtain a given lift. The distribution of induced velocity 
is, however, rather different from the calculated values, 
particularly in the central portion of the wake. In some 
cases, the calculated induced-velocity field in the usual 
location of short-span tail surfaces could be incorrect by 
as much as 1.5z? 0 . 

NATURE OF THE FLOW 

At this point the reason for the rather severe discrepancies 
between the calculated and the measured induced velocities 
in the rearward portion of the flow is considered. Refer- 
ence 1 has already noted the rolling up of the vortex sheet 
behind the rotor. Figures 39 to 41 examine this phenome- 
non more thoroughly for three different flight conditions. 
These figures are contour maps of the measured local mean 
vorticity behind the rotor. The intersection of the outer- 
most edge of the assumed wake (refs. 3, 4, and 6) with the 
survey plane is showui in each case by the dashed-line 
ellipse. If the assumed wake vortex, pattern is correct, 
the contours of vorticity should be along this wake edge 
(or within it, depending upon the disk-load distribution). 

In general, this is not true since the dominant feature in 
each case is the vorticity associated with the well-developed, 
rolled-up vortices behind the rotor tips. It should be noted 
that these vortices have only been displaced vertically 
about half as far as the center of the wake. The lower 
half of the elliptical sheet is essentially completely rolled 
up at the trailing edge (figs. 39, 40 (a), and 41). The upper 
half of the elliptical sheet (originating from the rear of the 
disk) is still visible at X/R=-~ 1.07, but it also is completely 
rolled up shortly after that point (figs. 40 (b) and 40 (c)). 

The similarity between the flow behind the rotor and 
that behind a low-aspect-ratio wing is very close and it is 
possible to draw a somewhat more quantitative analogy. 

Reference 15 examines the rolling up of the trailing-vortex 
system behind a low-aspect-ratio wing. For lift coefficients 
of the order of several tenths, the roll-up is essentially com- 
plete in a very short distance. Figure 9 of reference 15 
shows that, for an untwisted, rectangular wing of aspect 



Figure 39. — Distribution of vorticity behind rotor. X/R= 1.07; 
x=75.0°; m=0.095. 





6 1 2 10 8 6 4 2 0 -2 -4 -.6 -8 -1.0 -1.2 -14 




Advancing side Y/R Retreating side 


(a) X/R= 1.07. 

(b) X/R= 2.07. 

(c) X/R=3A4. 

Figure 40. — Distribution of vorticity behind rotor. x=82.3°; 

/u=0.140. 



Advancing side Y/R Retreating side 


Figure 41. — Distribution of vorticity behind rotor. X/R= 1.07; 
X = 83.9°; /*= 0.232. 
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ratio 4/7r (the same aspect ratio as a rotor), and for the lift 
coefficients used herein (0.720, 0.373, and 0.122), the roll-up 
distances should be 0.33, 0.64, and 1.97 times the chord 
length, respectively. It is more difficult to determine the 
rate of roll-up for something as complex as a rotor, but 
figures 39 to 41 indicate qualitatively that it is also extremely 
rapid. 

These observations are further substantiated by the 
induced-velocity measurements previously examined in 
figures 19 and 21. The measured data in these figures show 
a marked discontinuity (which would be expected at the 
edge of an elliptical wake) only in the vicinity of the upper 
wake edge at X/R=\.07. Farther behind the rotor 
(A/R= 2.07 and 3.14), there is no pronounced discontinuity; 
this condition indicates the absence of a strong vortex sheet 
in this region. 

Figures 22 (f) and 22 (g) seem to show a discontinuity 
somewhat above the upper edge of the assumed wake. 
The previous analogy with low-aspect-ratio wings indicates 
the possibility that the roll-up would not be as rapid in 
this case because of the low rotor lift coefficient (C L = 0.122). 
The difference between the calculated and the measured 
flow fields may be explained on this basis. In the forward 
portions of the flow, the elliptical wake has not rolled up, 
and it therefore resembles the wake-vortex model used in 
the calculations. The wake here, has, of course, the greatest 
effect on the induced velocities in the forward portions of 
the flow and essentially determines their magnitude. Con- 
sequently, in this part of the flow, calculations give a reason- 
able picture of the actual induced velocity. 

Farther rearward, however, the distorted portion of the 
wake has a greater effect in determining the actual induced 
velocity. The calculated flow field, therefore, is increasingly 
more inaccurate as its distance downstream increases. 

A short distance behind the rotor, it should be possible 
to assume that the flow will be the same as that for an 
equivalent wing. The minimum distance behind the rotor 
at which this procedure will be valid is a function of the 
lift coefficient. 

From reference 16 the induced velocity in the plane of 
symmetry of a uniformly loaded rectangular wing is 
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If it is assumed that the quarter chord of the equivalent 
wing and the lateral center line of the rotor coincide and 
also that s=0.S5R (the approximate center of gravity of 
the vortex system), 
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If nondimensionalized, equation (4) becomes 
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Now 


r= 


2 P VR 


and, for skew angles near 90°, 
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Substituting equation (6) into equation (5) and dividing 
both sides of the resulting equation by equation (7) gives 
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Figure 42. — Comparison in the plane of symmetry of measured 
values of induced-velocity ratio u/t’ 0 with those calculated for an 
equivalent wing. x=82.3°; 0.140. 
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The flow calculated by this equation is symmetrical about 
Z/R= 0 ; that is, it does not allow for the vertical displacement 
of the wake. As pointed out previously, the centers of the 
vortices are displaced downward an amount 

A Z 1 X . 

A R ~~ 2 R cot X 

The induced velocity for this assumed wing has been 
calculated from equations (8) and (9) and is compared with 
the flow measured at Ayi?=2.07 and 3.14 (for /x=0.140 and 
X=82.3°) in figure 42. As expected in these rearward 
planes, it is a closer approximation to the actual flow than 
are the flows calculated on the basis of a skewed cylindrical 
wake (fig. 21, for example). 

CONCLUSIONS FOR PART III 

The results of this wind-tunnel investigation of the in- 
duced flow near a single rotor in simulated cruising and high- 
speed flight are as follows: 

1. As far rearward as three-quarters of a diameter behind 
the leading edge of the rotor disk, the normal component 
of induced velocity of a lifting rotor may be calculated 
with good accuracy by available theory provided that a 
realistic nonuniform disk loading is assumed. 

2. Rearward of this three-quarter diameter point, cal- 
culation of the induced velocity becomes increasingly inaccu- 
rate because the trailing-vortex system rolls up very rapidly, 
and the wake vortex pattern behind the rotor is altered 
from that assumed as a basis for the calculations. 

3. At locations well behind the rotor, the induced flow 
may be calculated more accurately by assuming the rotor 
to be a uniformly loaded rectangular wing, provided that 
the rotor lift coefficient is sufficiently high. 

4. There is little effect of the assymmetries associated 
with tip-speed ratio on the induced flow along the longi- 
tudinal axis of the rotor. 

5. The induced-velocity calculations of Mangier and 
Squire can be modified to account for the wake crossing 
the flow because of the finite lift coefficient. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 7, 1956. 
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APPENDIX 

BASIC DATA 

At present the knowledge of the flow field near a lifting 
rotor is very limited. Because of this lack of information 
substantially all of the basic data obtained in the experi- 
mental investigation are presented in this appendix. In each 
case, the rotor flight condition is identified b}^ the skew angle 
X and the tip-speed ratio n. Table I (p. 27) identifies the 
flight conditions more completely. 

Figures 43 to 45 show the measured stream angles in the 
flow; figures 46 to 48 show contours of dynamic-pressure 
ratio; and figures 49 to 69 give the measured values of the 
induced velocity ratio. The fairings of the measured data, 
shown in figures 49 to 69, are those used in the analysis of 
part III. 


48 


REPORT 1319 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Deviation, deg 

(a) X/R=-1.0. 

Figure 43.— Measured stream angles in vicinity of rotor. Horizontal and vertical components of vectors represent yaw and pitch angles. Base of vector is at point of measure- 

ment. x=75.0°; ^=0.095. 
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Figure 43. — Concluded. 
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Figure 44. — Measured stream angles in vicinity of rotor. Horizontal and vertical components of vectors represent yaw and pitch angles. Base of vector is at point of 

measurement, x— 82.3° 0.140. 
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(f) A7tf = 2.07. 
Figure 44.— Continued. 
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(g) XIR = 3.14. 
Figure 44. — Concluded. 
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Figure 45. — Measured stream angles in vicinity of rotor. Horizontal and vertical components represent yaw and pitch angles. 

x = 83.9°; ju = 0.232. 
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(a) X/R=- 1.0. 

Figure 46. — Contours of dynamic pressure q/qo in the vicinity of the rotor. x = 75.0°; ^ = 0.095. 
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(c) X/R = 3.14. 
Figure 47. — Concluded. 
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(a) X/fl=1.07. 

(b) X/R = 2.07. 

(c) X/R = S. 14. 

Figure 48. — Contours of dynamic pressure q/q 0 in the vicinity of the rotor. x=83.9°; n = 0.232. 
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Figure 49.- — Measured values of induced-velocity ratio v/v 0 . X / /?= — 1.0; x = 75.0°; m — 0.095. 



INDUCED VELOCITIES NEAR A LIFTING ROTOR WITH NONUNIFORM DISK LOADING 


73 




Figure 52. — Measured values of induced-velocity ratio i>/t>o- X/R= 0.5; x = 75.0°; ^ = 0.095. 
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Figure 53. — Measured values of induced-velocity ratio v/v 0 . X/R= 1.07; x=75.0°; /u = 0.095. 
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Figure 54. — Measured values of induced-velocity ratio r/% X/ft =2.07; x=75-0°; ju = 0.095. 
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Figure 55. — Measured values of induced-velocity ratio v/v 0 . X/R= 3.14; x = 75.0°; ^=0.095. 
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Figure 56. — Measured values of induced-velocity ratio v/v 0 . X/R= — 1.0; x=82.3°; m=0.140. 



Advancing side Y/R Retreating side Advancing side Y/R Retreating side 


Figure 57. — Measured values of induced-velocity ratio t>/^o. X/R- — 0.5; x = 82.3°; m = 0.140. 



Figure 58. — Measured values of induced-velocity ratio v/vq. X/R= 0; x = 82.3°; /*=0.140. 
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Figure 59. — Measured values of induced-velocity ratio v/v 0 . X/R = 0.5 ; x=82.3°; m= 0.140. 
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Figure GO. — Measured values of induced-velocity ratio v/vq. X/R= 1.07; x=82.3°; m=0.140. 
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Figure 61. — Pleasured values of induced-velocity ratio v/v 0 . X/ R = 2.07 ; x=82.3°; ^ = 0.140. 
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Figure 62. — Measured values of induced-velocity ratio v/v 0 . X/R = 3.14; x=82.3°; /* = (). 140. 
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Figure 63. — Measured values of induced-velocity ratio v/v 0 . X/R= — 1.0; x=83.9°; 

m = 0.232. 
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Figure 64. — Measured values of induced-velocity ratio v/v 0 . X/R= — 0.5; x = 83.9°; 

m = 0.232. 
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Figure 65. — Measured values of induced- velocity ratio v/v 0 . X/R = 0; x = 83.9°; 

u= 0.232. 
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Figure 66. — Measured values of induced- velocity ratio v/vo. X/R= 0.5: 

X=83.9°; ju = 0.232. 


86 


REPORT 1319 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 




Figure 67. — Measured values of induced-velocity ratio v/vo. X/R= 1.07; x=83.9°; n=0.232. 
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Figure 69. — Measured values of induced-velocity ratio vj v 0 . XJR = 3.14; * = 83.9°; ^ = 0.232. 
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